Interference of Topologically Protected Edge States 
in Silicene Nanoribbons 



Motohiko Ezawa^ and Naoto Nagaosa^'^ 
^Department of Applied Physics, University of Tokyo, Hongo 7-3-1, 113-8656, Japan and 
^ Cross Correlated Materials Research Group (CMRG) and Correlated Electron Research Group (CERG), 

ASI, RIKEN, Wako 351-0198, Japan 

Silicene is a two-dimensional quantum spin-Hall insulator. We study the edge channels of silicene nanorib- 
bons from the viewpoint of the topological protection and the interference between the two edges. It is found 
that the behaviors of the helical edge channels (HECs) are completely different between the armchair and zigzag 
edges. The penetration depths ^ of the HEC is antiproportional to the spin-orbit gap for the armchair edge 
^arm ^ Hvf / A (vf'. thc Fcrmi velocity, A: the gap due to the spin-orbit interaction), while it remains as short 
as the lattice constant for the zigzag edge. Zero-energy states disappear in armchair nanoribbons due to an 
interference of two edge states, while they remains in zigzag nanoribbons even if the width W is quite narrow. 
The gap 6 of HECs behaves as ^ hv^/W for A = 0, and as ^ ^ A exp[— VK/^arm] for A / for armchair 
nanoribbons, while it remains essentially zero irrespectively of A for zigzag nanoribbons. 



Graphene, a monolayer honeycomb structure of carbon 
atoms, is one of the most fascinating topics in condensed 
matter physics^. Graphene has two Dirac fermions at the K 
and points in the Brillouine zone, which govern the low 
energy properties and produce many nontrivial phenomena. 
Another intriguing aspect of graphene is the edge channels. 
There are two types of edges2rd, i.e., zigzag and armchair 
edge [Figs. 1(a) and (b)]. The edge channel is strictly local- 
ized at the outermost atoms for the zigzag edge, while there 
is no exponentially localized channel for the armchair edge. 
For the zigzag edge, the flat dispersion appears connecting the 
two valley K and points [Fig.2(al)], while only one val- 
ley is relevant to the low energy states for the armchair edge 
[Fig.2(bl)]. In the armchair edge case, the finite gap S appears 
due to the finite- width effects, i.e., S ~ v^/W with the Fermi 
velocity vip and the width W. For the zigzag edge, the gap is 
always zero for the edge channel even for quite narrow width 
W since it is exactly localized at the outer most atoms. 

We have up to now neglected the spin-orbit interaction 
(SOI). With the SOI, graphene turns into a quantum spin-Hall 
(QSH) insulator^, although the QSH effect can occur only at 
unrealistically low temperature^'^ due to its too small SOI. Re- 
cently, it has been shown^ that the QSH effect is naturally 
realized in a honeycomb structure of silicon named silicene, 
which was experimentally manufactured^ last year. Sil- 
icene has enormously rich physics ^^^^ in view of topologi- 
cal insulators. Therefore, silicene offers an ideal laboratory to 
study the topological properties of honeycomb lattice system. 

According to the bulk-edge correspondence there 
should be gapless helical edge channels when the bulk 
states are topologically nontrivial. Therefore, once the SOI 
is introduced, the edge channels of graphene without the SOI 
turn into the topologically protected helical edges. It is an im- 
portant issue to reveal how this crossover occurs as the SOI 
increases, which we address in the present paper. We found 
the SOI and the topology play completely different roles be- 
tween the armchair and zigzag nanoribbons, though the topo- 
logically protected helical edge channels appear in both cases. 
The penetration depth of the edge channels and the conse- 
quent hybridization^^'^^ between the two edges are examined. 
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FIG. 1: (Color online) Illustration of (a) zigzag and (b) armchair 
nanoribbons. Its width W is defined by the number of hexagons in 
a unit cell. Here we have taken = 5 for zigzag and = 8 for 
armchair nanoribbons. (c) The hexagonal Brillouin zone. The states 
near the Fermi energy are tt orbitals residing near the K and K' 
points at opposite corners of the hexagonal Brillouin zone. The bulk 
band structure of nanoribbons are obtained by projecting the band 
structure of the bulk from the direction depicted in the figure. The K 
and points are identified in the armchair edge. 



and they remain essentially zero for the zigzag nanoribbons in 
sharp contrast to the armchair one. 

Hamiltonian: We employ the Kane-Mele model^ to de- 
scribe silicene as given by 
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^a/3Cj/3, (1) 



where c]^ creates an electron with spin polarization a at site 
z in a honeycomb lattice, and (z, j) / {{i^j)) run over all the 
nearest/next-nearest-neighbor hopping sites. The first term 
represents the usual nearest-neighbor hopping with the trans- 
fer energy t, while the second term represents the effective 
SOI with Aso, where a = (da^, cr^, cr^) is the Pauli matrix of 
spin, with = +1 if the next-nearest-neighboring hopping 
is anticlockwise and Uij = — 1 if it is clockwise with respect 
to the positive z axis. 

Tthe Hamiltonian ([T]) describes the basic nature of sil- 
icene, that is a honeycomb structure of silicon atoms, where 
t = 1.6eV and Aso = 3.9meV^'2^ The SOI Aso is reason- 
ably large. The crucial advantage enjoyed by silicene is its 




FIG. 2: (Color online) Band structure of (al,a2,a3) zigzag and 
(bl,b2,b3) armchair nanoribbons. We have taken Aso = for 
(al,bl), Xso/t = 0.1 for (a2,b2), and Xso/t = 0.2 for (a3,b3). We 
have taken the width W = 16. The vertical axis is the energy in 
unit of t and the horizontal axis is the momentum k. Zero-energy 
edge modes are present in the zigzag nanoribbon (a2,a3) but not in 
the armchair nanoribbon (b2,b3), though the bulk is a topological 
insulator for Aso / in both cases. The cyan region (red curve) 
represents the band of the bulk (edge). The bulk spectrum takes the 
minimum at the K and K' points. (The K and K' points are identified 
in the armchair nanoribbon.) The bulk mode is well described by the 
analytic formula (|5]), while the edge mode by the dispersion relation 
iT2\ for a zigzag nanoribbon and ([T4t for an armchair nanoribbon. 
The band gap of the bulk (edge) is denoted by 2A (2^), which in- 
creases (decreases) as Aso increases. 



buckled structure separating the sublattice planes for A sites 
and B sites by a distance 2£ — 0.46 A. It generates a staggered 
sublattice potential oc 2iEz between silicon atoms at A sites 
and B sites in electric field E^!^^ . Furthermore, we may gen- 
erate the Haldane interaction^^ term with strength A^ by way 
of photo-irradiation^^. It is also possible to include the stag- 
gered exchange magnetization^^ with strength AM. They are 
summarized as an additional term AJY to the Hamiltonian ([T]), 



+AMy^/ii4^(7^ 



(2) 



where /i^ = ±1 for i representing the A (B) site. This ad- 
ditional term provides silicene with enormously rich physics. 
In the present problem we are able to control the penetration 
depth ^ of the edge mode experimentally by tuning these ex- 
ternal parameters. 
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FIG. 3: (Color online) Band gap 28 of armchair nanoribbons as a 
function of the width W for various spin-orbit interactions Aso jt — 
0, 0.05, 0.1, 0.15, 0.2 (from up to down). The vertical axis is the 
energy in unit of t and the horizontal axis is W . 



For the sake of simplicity we numerically investigate 
nanoribbons based on the Kane-Mele model ([T]), but we in- 
clude the term (|2]) in constructing the low-energy theory based 
upon which we carry out an analytic study of nanoribbons. 
We define the width W of the nanoribbon as the number of 
hexagons in a unit cell as shown in Fig.l. The unit cell con- 
tains 2iy + 4 (2iy + 2) silicon atoms for armchair (zigzag) 
nanoribbons. We have diagonalized numerically the Hamil- 
tonian ([T]) to obtain the eigenvalues and the eigenstates, from 
which we find the band structure and the wave function. 

In FiglJl we show the evolution of the dispersions of the 
electronic states as the SOI is increased for zigzag [(al)-(a3)] 
and armchair [(bl)-(b3)] nanoribbons, respectively. There are 
two types of gaps, one for the edge part (2^) and the other for 
the bulk part (2 A). In the case of zigzag nanoribbon, the flat 
dispersion begin to cant and yield the helical modes intrinsic 
to the QSH insulator as the SOI is introduced [Figl2ta2,a3)]. 
It is noted that the gap 28 for the edge channel does not ap- 
pear even with finite SOI although the gap 2 A becomes finite 
for bulk states. On the other hand, in the case of armchair 
nanoribbons, 28 due to the finite size effect decreases as the 
SOI is increased. [FigI3bl)-(b3)]. Below, we analyze the 
wavefunctions of the edge channels in more details to examine 
the different behaviors of zigzag and armchair nanoribbons. 

Band structure of Silicene nanoribbons: We show the 
width dependence of the band gap 28 for several fixed values 
of Aso in FiglJl The band gap oscillates in the period of three, 
as is a well-known feature^ of armchair nanoribbons. When 
Aso = (Aso 7^ 0), the band gap decreases antipropotionally 
(exponentially) as W increases for ModsW^ ^ 0. 

We show the absolute value of the real-space wave function 
in FiglH When Aso = 0, the wave function is constant for 
ModsVK = and almost constant for ModsVK ^ across 
the nanoribbon. The peaks emerge at the both edges as Aso 
increases. They are the zero-energy edge modes required by 
the bulk-edge correspondence, as we shall soon demonstrate 
based on analytic formulas. We remark that there is a consid- 
erable amount of overlap between them. The overlap becomes 
smaller as Aso increases. The order of the overlap is measured 
by the penetration depth ^ of the edge mode. 

Low-energy Dirac theory: We proceed to construct the 
low-energy theory to make a further study of the zero-energy 
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FIG. 4: (Color online) Real- space wave functions of armchair 
nanoribbons for SOI Xso/t = 0,0.05,0.1,0.15,0.2. We take (a) 
W = 48 and (b) W 49. Note that mod348 and mod349 1. 
The vertical axis is the energy in unit of t and the horizontal axis is 
the y-axis of the nanoribbon with the width 2L. 



FIG. 5: (Color online) The bonding state ^ obtained analytically 
presents a good fit of the numerically determined wave function. 
Here we take examples of = 48 and 49. The vertical axis is 
the energy in unit of t and the horizontal axis is the y-3xis of the 
nanoribbon with the width 2L. 



modes and their overlap in a nanoribbon. We adopt the Hamil- 
tonian H + AH in order to apply our results to a realistic ma- 
terial such as silicene. The low-energy theory in the Krj (K or 
K^) valley is given by the Dirac Hamiltonian, 



Hrj = hvv (rikxTx + kyTy) + XsoCTzV'^z 



(3) 



where v^^ = ^at = 5.5 x lO^m/s is the Fermi velocity with 
the lattice constant a = 3. 86 A, and is the Pauli matrix of 
the sublattice pseudospin. This Hamiltonian describes a four- 
component Dirac fermion indexed by the spin = ±1 and 
the pseudospin = ±1 for each valley r] = ±1. 

The coefficient of is the mass of Dirac fermions in the 
Hamiltonian ([3]), 



(4) 



which plays the most important role in the physics of sil- 
icene. It is intriguing that the spin- valley dependent mass A^^ 
may be positive, negative or zero. A nontrivial topological 
charge is generated when A^^ takes a negative value. Silicene 
is shown to be a QSH insulator without the external fields 
(E^ = 0, = 0, AM = 0). The band gap is given by 
2 1 A^ I . The energy spectrum reads 



E{k) = ±J{hv^)^k^^{Al)^, 



(5) 



which is illustrated by taking A^^ = Aso in Fig 121 It gives a 
good approximation to the band structure of the bulk. 

Low-energy theory of armchair nanoribbons: We inves- 
tigate the zero-energy edge modes of armchair nanoribbons. 
We take the x direction as the translational direction of a 
nanoribbon. The zero-energy edge modes appear at k^ = 0. 
The transverse momentum ky is determined by solving the 
zero-energy solution of (0), £^(0) = with (|5]). It is solved as 

ky = ±i\Al\/ihv^). (6) 

The wave function for the edge located at ±L reads 



±- 



ivTL) 



(7) 



up to a normalization constant, where Q{\y\ — L) = 1 for 
\y\ < L and (d{\y\ - L) = for \y\ > L. Here, L and W are 
related as L = ^Wa. The penetration depth is given by 



hv,/\Al 



(8) 



We have demonstrated the emergence of the zero-energy 
modes © at the two edges (y = ±L). Their wave functions 
mix due to the interedge interaction, and form the bonding 
state as the ground state. The bonding state is given by 



(9) 



up to a normalization constant. The wave functions © and (|9]) 
present remarkably good approxiamtions to the envelop func- 
tions of the numerically calculated wave functions [Fig 121. 

The bonding state is no longer a zero-energy state due to 
the mixing. The energy is estimated as 
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r-L{y)^+L{y)dy = |A^Jexp(-2L/ea. 



(10) 

as the overlap of the two edge states. This is the reason why 
zero-energy edge modes disappear from the energy spectrum 
of armchair nanoribbons. 

The effective Hamiltonian of the armchair edge states reads 



H = a,Tf^'hv^k,^STf^', 



(11) 



where r^^^^ is the Pauli matrix for the edge pseudospin, 
Tz^^ = ±1 for the top and bottom edges. The first term de- 
scribes the two edge states (r^^^^ = ±1) with the opposite ve- 
locity each of which carries the up and down spins (az = ±1). 
The second term describes the mixing of the two edge states. 
By diagonalizing the Hamiltonian ([TT]) . the eigenvalue is 



E = ±^{hvF)^kl^S^. 



(12) 



This gives a good approximation of the edge mode in Fig 121 
where S ^ S. Namely, the overlap integral produces the 
gap of the edge states. We conclude that, strictly speaking, 
zero-energy edge modes never appear in armchair nanorib- 
bons. However, practically they appear provided the condition 
^arm <C L is Satisfied, where S 0. 
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FIG. 6: (Color online) Real- space wave functions of zigzag nanorib- 
bons for SOI Xso/t = 0,0.05,0.1,0.15,0.2. They are well de- 
scribed by the analytic formula ([TSt . Inset: A logarithm plot of 
the wave functions. Clearly it decreases linearly as the position in- 
creases. We take W — 24. 



Low-energy theory of zigzag nanoribbons: Finally, we 
investigate the wave functions of zigzag nanoribbons with 
the SOL We show the absolute value of the real-space wave 
function in Fig IS] The edge state in one edge is completely 
localized at A sites while the other edge at B sites when 
Aso = 0. The two states localized at the two edges are orthog- 
onal to each other, and there is exactly no overlap between the 
two edge modes. Furthermore, the penetration depth is zero. 
When Aso 7^ 0, although the totally localized state is not an 
exact solution, the wave function is almost localized at the 
edge. The overlap between the two edge states is found to be 
zero as well within the accuracy of our numerical calculation. 

The edge mode crosses the Fermi energy kx = tt as in 
FigI2a2,a3). Since the Dirac Hamiltonian © describes solely 
the low-energy theory near the K and K' points separately, 
it does not provide us with the low-energy Hamiltinian of a 
zigzag nanoribbon connecting the tips of the two Dirac cones. 
Nevertheless we are able to write down the phenomenological 



Hamiltonian for the zigzag edge states, 

by requiring a linear dispersion, 

E = ±\sohvYkx/t, 



(13) 



(14) 



as is the result of numerical analysis. The electron velocity in 
the edge states is almost constant and proportional to the SOI 
Aso- This dispersion gives an excellent fitting of the zero- 
energy edge mode as in FigI2a2,a3). 

The wave function of the zero-energy state is well fitted by 

^±L{y) = e{\y\ - L) exp [±{y T L)/U^] , (15) 

where the penetration depth is approximately given by 

Czig^alA^J/t. (16) 

It is interesting that the penetration depth of zigzag and arm- 
chair edges have an inverse relation, ^arm^zig — ahv^/t. 

Discussion: We have explored numerically and analytically 
the properties of the edge mode in silicene nanoribbons. The 
gap of armchair nanoribbon with the width of the order of 
10/im is lOmeV, which can be observed experimentally. We 
have found that the zero-energy edge modes is robust against 
the hybridization between the two edges in zigzag nanorib- 
bon, while the finite gap 5 in the edge channel is observable 
in the armchair nanoribbon. This offers an interesting possi- 
bility to construct the ideal situation where only the Coulomb 
interaction is effective between the two helical edge channels 
without the hybridization, where a new electronic liquid state 
has been proposed^^. 
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